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Abstract 

We consider a system of weakly coupled singularly perturbed semilinear elliptic 
equations. First, we obtain a Lipschitz regularity result for the associated ground 
energy function E as well as representation formulas for the left and the right 
derivatives. Then, we show that the concentration points of the solutions locate 
close to the critical points of E in the sense of subdifferential calculus. 
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1 Introduction and main results 

In the asymptotic analysis of the singularly perturbed elliptic equation 

- e 2 Au + u = f(x,u) inE", u>0 in R™, (P e ) 

there are well known situations where the associated ground energy function £ 
(cf. [30]) is C 1 -smooth and around its nondegenerate critical points the solutions 
u e of P e exhibit a spike-like profile as e goes to zero. This is the case, for instance, 
for the power nonlinearity 

f{x,u)=K(x)u«, K?<g, n>3, 

where K{x) is a suitable C 1 function (see e.g. [3,4] and references therein). It turns 
out that the C 1 (and higher) smoothness of S is related to the crucial fact that, for 
every fixed z 6 R™, the limiting autonomous equation 

- Ait + u = f(z, u) inR", m>0 in R", (P ) 

admits a unique solution, up to translations [30]. However, unfortunately, the 
uniqueness feature for Pq is a delicate matter and it is currently available only un- 
der rather restrictive assumptions on / (cf. e.g. [16]). What it is know, in general, 
is only that E is a locally Lipschitz continuous function which admits representa- 
tion formulas for the left and right derivatives (cf. [30, Lemma 2.3]). Motivated by 
these facts, recently, some conditions for locating the concentration points for P e 
in presence of a more general nonlinearity /, not necessarily of power type, have 
been investigated (see [23] and also [24]). The underlying philosophy is that when 
the limit problem P Q lacks of uniqueness up to translations, then the ground energy 
function S could loose its additional regularity properties. 

Nevertheless, in this (possibly nonsmooth) framework, it turns out that a nec- 
essary condition for the solutions u e to concentrate (in a suitable sense) around a 
given point z is that it is critical for E in the sense of the Clarke subdifferential dc , 
that is G dc^{z), or in a even weaker sense. The main theme of this note is the 
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search of suitable conditions for locating the spikes, as e — > 0, of the solutions to 
the semilinear model system 

'-e 2 Au + u = K(x)v q , inR", 

< -e 2 Av + v = Q(x)vP, inR", (S e ) 
u,v>0, inR", 

where p > 1 and q > 1 are lying below the so called "critical hyperbola" 

tfn = {fog) € (l.oo) x (l.oo): ^ + ^ = 1 - \ }, n > 3, 

which naturally arises in the study of this problem and constitutes the borderline 
between existence and nonexistence results (cf. e.g. [10,15]). 

Now, according to what discussed right above, the interest in looking for con- 
ditions for the spike location of the solutions to \S E \ is mainly motivated by the 
following simple observation: contrary to the scalar case, there is no uniqueness 
result available in the literature for the (radial) solutions to the (limiting) system 
associated with \S F l 

-Au + u = K (z)v q , inR", 

-Av + v = Q(z)vP, inR", {S z ) 
u,v>0, inR", 

where z 6 R" is frozen and acts as a parameter. As a consequence, in the vectorial 
case, we do not know whether the (suitably defined) ground energy map £ associated 
with \S e \ (cf. Definition 11.21) is C x -smooth and admits an explicit representation 
formula. Hence, the necessary conditions in terms of Clarke subdifferential (or 
weaker) appear here even more natural than in the case of a single equation. See 
Section 11.21 for the statements of the main results, Theorems 11.11 and 11.21 As far 
as we are aware, other criteria for the concentration have been established so far, 
but all of them consider the scalar case. We refer the reader e.g. to [3,20,29,30] 
for the case of power-like nonlinearities and to [23, 24] for more general classes of 
nonlinearities. 

Semilinear systems like jS E l naturally arise in the study of various kinds of 
nonlinear phenomena such as population evolution, pattern formation, chemical 
reaction, etc., being u and v the concentrations of different species in the pro- 
cess (see also [32] and references therein). Visibly, the interest in the study of 
the various qualitative properties of \S £ \ has steadily increased in recent times. 
In a smooth bounded domain f2, (Si) was extensively studied by Clement, Costa, 
De Figueiredo, Felmer, Hulshof, Magalhaes, van der Vorst in [10-12,14,15]. The 
asymptotic analysis with respect to s has been very recently performed both with 
Dirichlet and Neumann boundary conditions by Pistoia-Ramos [18,19] and Ramos- 
Yang [22]. In the whole space R ra , the existence of least energy solutions to jS e [ 
has been investigated by Alves-Carriao-Miyagaki, De Figueiredo, Yang and Sirakov 
in [1,2,13,26,27,31], whereas the asymptotic behavior with respect to e has been 
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pursued by Alves-Soares-Yang in [2] . Finally, for the exponential decay, the radial 
symmetry and the regularity properties of the solutions to \S z l , we refer the reader 
to the quite recent achievements of Busca-Sirakov and Sirakov [6,27]. 

The outline of the paper is as follows: in Sections ll.lll.2l we provide preliminary 
stuff such as the (dual) variational framework and the (dual) ground energy function 
S and we state the main results of the paper. Throughout Section [21 we deal with 
the Lip loc regularity and the representation formulas of the directional derivatives 
for S. Finally, in Section |3 we end up the proofs of the main results. 



1.1 The dual variational functional 

As it is known, if e.g. p and q are both less than ^r§, then system \S E l admits a 
natural variational structure (of Hamiltonian type) which is based on the strongly 
indefinite functional f e : H 1 (W n ) x iy x (R") -> R, 



/*(«,«)=/ s 2 Wu-Vv + uv~-^ f K{x)\v\ q+1 - — pj- f Q(x) 

JR™ it" JK" 



u\ p+1 . 



However, as already done in [1,2], for our purposes, as well as for dealing with 
possibly supercritical values of p or q, we consider a corresponding dual variational 
structure, mainly relying on the Legendre-Fenchel transformation (see e.g. [8,9,17] 
and references therein). In the following, we just briefly recall some of the core 
ingredients, referring to [1, Section 2] for expanded details on this framework. For 
^q-j- + > consider the linear operators 

T x : L^fR") W 2 ' 3 ^r{W l ) ^ L p+1 {R n ), 
T 2 : L^p 1 (R n ) -> W 2 ^{R n ) <-+ L q+1 {R n ), 

defined as 

Ti =T 2 = (-A + Id)- 1 . 

Notice that T\ and T 2 are continuous. Then, we consider the linear operator (take 
into account the proper Sobolev embeddings) 

T : L^iW 1 ) x L^r{R n ) -> L p+1 {R n ) x L q+1 (R n ), T = 



Ti 
T 2 



explicitly defined by 

( T v-,0 = i\Tim + &T 2 m, v?7 = (771,772), Vf = (£i,6). 

Finally we introduce the Banach space (Jf?, \\ ■ ||jr), 

L P (K") L 1 (R™) 

and the (dual) C 1 functional J e : — > R defined as 

P 1 ml p q f \m\ q 1 
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If if = (??i,?7|) is a critical point of J E , then (it e (x), v e {x)) = (u e (|), u e (f )), with 

(u e ,i; e ) = (Tnjf,T 2 »£) £ W 2 ^nL p+1 x lf 2 'T i nL ?+1 , (1.1) 

corresponds to a solution to ]S e j with u e (x), v e (x) — ► for |sc| — ► oo (see [1, p.677]). 
In light of the above summability, we have f e (u s , v e ) £ R for all e > 0. Analogously, 
associated with d5* z | ) , we introduce the limiting functional I 2 : Jf° — > K 

From the viewpoint of our investigation, the main advantage of exploiting the 
dual variational functional I z is that it admits a mountain-pass geometry and the 
mountain-pass value corresponds to the least possible energy of system \S z l - As 
we shall see in the next section, this allows to provide in the vectorial framework 
a suitable definition of ground energy function with nice features, similar to those 
available in the scalar case. 

1.2 Preliminaries and the main results 

In order to state the main achievements of the paper, we need some preparatory 
stuff. For the sake of self-containedness we shall also recall a few pretty well known 
notions from nonsmooth calculus (see e.g. [7]). 

Definition 1.1 Let / : R™ — > R be a locally Lipschitz function near a point z £ R n . 
The Clarke subdifferential of / at z is defined by 

d c f{z) := {q £ R" : f(z, w)>r]-w, for every w £ R' 1 }, 

where f°(z,w) is the generalized derivative of / at z along w £ R'\ defined by 

j [z; w) :— nmsup . 

Definition 1.2 The (dual) ground energy function £ : R n — > R of \S z l is given by 

£(*):= inf. I z (rf), 

where J\f z is the Nehari manifold of I z , that is 

K = {V S J? : 7? 7^ (0, 0) and I' z ( V )[n] = 0}. 
We shall denote by /C C 1" the set of Clarke critical points of E, namely 

K := {z £ R™ : £ d c Z{z)}. 
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Definition 1.3 We say that the pair (u e , v E ) is a strong solution to system jS e [ if 
it is a distributional solution and (u e ,v s ) G W 2,( - q+1 ^ q (R n ) x W^ 2 '(P +1 )/P(M n ). We 
say that the pair if — (rjf, ryf) corresponding to (u E , v E ) through is the related 
dual solution. 

Definition 1.4 We set 

£ := {z G M™ : there exists a sequence of strong solutions (u Eh , v Eh ) of \S E \ with 
\u Sh (z)\, \v Eh (z)\ > S for some S > 0, \u Eh (z + e h x)\, \v Eh (z + e h x)\ — > 
as |x| — * oo uniformly w.r.t. ft,, and Sh~ n f Eh (u Eh , v Eh ) — > E(z) as ft. — > oo}. 

We say that £ is the energy concentration set for \S E \ . 

Assume that K, Q G C^M 71 ) and 

a < A(x) < /3, a<Q(x)<f3, for all iel", (1.2) 

iV-K^a;)!, |VQ(a;)| < Ce M|x| , for all ig!" with |x| large. (1.3) 

for some positive constants a,f3,C and M. 

The main result of the paper, linking the energy concentration set £ with the 
set K. of Clarke critical set of E, is provided by the following 

Theorem 1.1 Assume that K,Q e C^K" ) and that lfO |l -lfO | ftofti Tften £ C }C. 

Remark 1.1 By [2, Theorem 1], under suitable assumptions, if there exists an 
absolute minimum (or maximum) point z* for E, then z* G £ ^ 0. 

Remark 1.2 As a straightforward combination of Theorem ll . ll with the well known 
convex hull characterization of dcE(z), if z is a concentration point for jS E [ , then 

G Co{ limVE(^) : £ ti and ^ -> z}, 

where Co{A} denotes the convex hull of A and fl is any null set containing the set 
of points at which E fails to be differentiable. 

Corollary 1.1 Under the (unproved) assumption that, for all z G R™ , system \S z l 
admits a unique positive solution (up to translations), E is C 1 -smooth and 

, q+l P+l . 

£ C Crit(Q~ Kp^), 

where Ciit(f) denotes the set of (classical) critical points of f. 

In the following definition we consider solutions which concentrate close to a 
point z, with bounded energy but not necessary stabilizing towards E(z). 
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Definition 1.5 Let m > 1. We set 

S m := {z £ R™ : there exists a sequence of strong solutions (u eh ,v eh ) of \S E l with 
\u Sh (z)\, \v £h (z)\ > S for some S > 0, |u £(1 (z + |v eh (z + £^ar)| -> 

as \x\ — > oo uniformly w.r.t. h, and e/j — n f Sh (u £h ,v Sh ) — > to as ft. — > oo}. 

We say that £ m is the concentration set for ]6' E | | at the energy level m. 

Definition 1.6 Let m > 1 and z S R n . For every io £ R™ we define rj m (w) by 



r z ,mH := sup 
tfeG m (z) 



' »76G m (z) 



p + ldw 7 R "QV(z) q+ldw J K " ifV( z ) 



where G m (z) denotes the set of all the nontrivial, radial, exponentially decaying 
solutions of jSzl having energy equal to m. 

It is readily seen that Tj m (w) £ R for all z, w in 1" (see the proof of (|2 . 1 If) ^ . It 
is also straightforward to check that, for any z £ R", the functions {w i— > rj m (w)} 
are convex. 

Definition 1.7 Let m > 1. We set 

Km := |z e r™ : o e ar- m (o) n ar+ m (o)}, 

where d stands for the subdiffcrcntial of convex functions, 

dTl m (0) = 6 R" : TT m H > £ • itf, for every w 6 R™ j. 

It is known by standard convex analysis that dTJ m (0) ^ 0, for every z £ R n . 
Observe that z £ K m if and only if is a critical point for both T~ m and T+ m . Of 
course, if G m (z) = {770} was a singleton, then z £ K, m if and only if 

dl 

r- m M = r+„» = ^(770) = 0, Vw £ R". 

Without forcing the energy levels of the solutions to approach the least energy 
of the limit system, we get the following correlation between the sets £ m and IC m . 

Theorem 1.2 Assume that K,Q £ C x (M n ) and JOJ-JOJ /iota. T/ien £ m C /C m . 



2 Properties of the ground energy function 

Before coming to the proof of the results, we need some preliminary stuff. 
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2.1 Some preparatory lemmas 

The next proposition is well known (see e.g. [31]); on the other hand, for the sake 
of completeness and self-containedness, we report a brief proof. 

Proposition 2.1 Let z <E E". Then (u,v) e W 2 ^{R n ) x W 2,S ¥~(R n ) is a solu- 
tion to ]S Z [ if and only if r\ — (771,772) = (T 2 ' 1 v,T^ 1 u) E Jf? is a critical point of 
I z . Moreover, there holds f z (u,v) — / z (771 , 772 ) ; where f z is the functional defined as 

/,(«,«)=/ Vu-Vv + uv--^ [ K{z)\v\" +1 - [ Q(z)\u\p +1 . 

Jr™ Jr™ Jr" 

Proof. Observe first that, if (u,v) G W 2 '*^ (W 1 ) x W 2 '^ {R n ) solves RQ , taking 
into account the Sobolev embedding, the value f z (u,v) is indeed finite (cf. I|l.lfl l. 
Let (u,v) be a solution to \S z l - Then, since 

m = t^v, m = ?V V 

we have 

77 2 = T^u = -Au + u = K{z)v«, 
7/1 = T~ x v = -Atj + v = Q{z) u p. 
Therefore, we get 

T 2 r]i = v = — ^ — and T1772 = u = — ^ — , (2.1) 

Ki{z) Qp(z) 

and so (771, 772) is a critical point of I z . Vice versa, if (771, 772) is a critical point of J 2 , 
it is readily seen that (|2.1[1 hold, so that (T1772, T2771) = (u,v) is a solution to ]S z j 
(cf. [1, p. 677]). Furthermore, on the solutions to jS z j , we have 

/»(«.«) = (3 " / ^K +1 + (2-^TT)/ K{z)v«+\ 

JR" JR™ 

Then, in light of 1)2. l[l . we have 

= (^t - 1) £ m T im + - 1) £ 772 r 2?7l 

= (^-1) / r „(- a - + ^+(^t-2) Jj- Au + U > 

= (|-^t)/ *K +1 + (l-^)/ iW +1 = /*(^), 
JK" Jr« 

which concludes the proof. 
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Definition 2.1 We say that r\ G is a dual solution to \S z l if it is a critical point 
of I z . We say that 77 is a dual least energy solution to jS z j if it is a dual solution 
and, in addition, I z (rj) = S(z). 

The next property, classical in the scalar case, will be pretty useful for our 
purposes. 

Lemma 2.1 For every z G R™, let us set 
b\(z) := inf sup / z (to), 

b 2 {z) := inf I z {n) = 

bs(z) := inf {l z (ri) : 77 G iff \ {0} is a dual solution to \S z l }. 

Then b\(z) = 62(2) = 63(2). Moreover {z £(2)} is continuous. 

Proof. The first equality follows from [1, Lemma 2]. Moreover in [1] it is proved 
that bi(z) — b2(z) is a critical value so that also &2(z) = b^(z) follows. Finally, by 
virtue of [2, Lemma 1], we know that £ is continuous. 

Lemma 2.2 Let z G W 1 and define the (nonempty) set 

JT+ := {r, e JT : / (T»j, »/)>()}. 

Then, for every n G J4?+, there exists a unique maximum point t n > of the map 
<j>: t G (0, oo) t — > I z (trj). In particular, t^n G M z . 

Proof. Let us observe that if 4>'(t) = 0, then 

(Tr l ,r,) = t^ [ ^+t^ ' ^ 



: {z) JR" Ki(z) 



i-p 



Since the function g(t) = At~~S~ + Bt~p~ with A, i? > is strictly decreasing for 
t > 0, then has at most one critical value. It is easy to see that for all 77 G Jrf?, 
4>(t) > for t small, while if r\ G it is readily seen that (f>(t) < for big /j's. 

2.2 Conjecturing the representation of £ 

Consider for a moment the equation 

-e 2 Au + V{x)u = K(x)u p , in E" , (2.2) 

with p subcritical and V and K potentials functions bounded away from zero. By 
the results of [16], we know that there is uniqueness (up to translation) of positive 
solutions for 

-Au + u = u p , in IT, 
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and, by a suitable change of variable, also for the "limit" problem at x = z of (|2.2[l 

-Ait + V(z)u = K(z)u p , inM n . 

This allows to give an explicit representation for the ground state function associated 
with Q2.2I) . merely depending on the potentials V and K (see for example [4,30]): 

p+l _ n 

y p-i 2/ z ) 

Z(z)=r \ 2 W , (2.3) 
Kp- 1 (z) 

for a suitable positive constant T. On the contrary, as already observed, up to our 
knowledge there is no (known) uniqueness result for the elliptic system 

-A£ + £ = C 9 , -A( + C-e p , inM", (2.4) 

and so, in general, we cannot provide an explicit expression for £ for \S z j . Slightly 
more in general, if V is smooth and a < V(x) < /3, consider the system 

-Au + V(z)u = K(z)v 9 , inM", 

-Au + V(z)v = Q(z)uP, inM™, (2.5) 
u, v > 0, in R". 

Assuming for a moment that 12.41) has a unique solution (£, Q, then we claim that 

(p+l)(q+l) n , „ 
V p,-l 2 ( Z ) 

sW = r— t±]~~> (2-6) 

Qp?- 1 (z) iif ps-i (z) 
for a suitable positive constant I\ Indeed, by rescaling 

u(x) = zui^(^lx) and v(x) = t^CC/-^); 



where we have set 



fi = n{z) := Vi(z), 

g+l 



GJ! = W\{z) 
-CD 2 = ZU 2 (z) 



Qpi-i [z) Kpi- 1 (z) ' 
p+l 
Vpi-i (z) 



it is easy to see that (it, u) is the unique solution of the system l|2.5|) . Hence, by 
a straightforward calculation, we reach l|2.6(l . Let us observe that the exponent of 
V{z) in (|2.6|) is equal to zero if, and only if, the pair (p, q) belongs to Then, 
for problems with powers p, q close to the set ^ , the potential V is expected to 
have a weak influence in the location of concentration points. Notice that the same 
phenomenon appears in the scalar case (cf. formula l|2.3|l ). since 2±_ — 2 ~ if and 
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only if p ~ 2±| = 2* — 1, where 2* is the critical Sobolev exponent for iJ 1 . Finally, 
we just wish to mention that, incidentally, the exponents 

p+1 „ q+1 



pq-1 pq-i 

in formula i|2.6|) also arise in the study of the blow-up rates for the parabolic system 

u t = Au + v q , v t =Av + u p , xEfl,t>0, 

with initial data u(x,0) = uq(x) > 0, v(x,0) = vq(x) > and Dirichlet boundary 
conditions u = v — on dfl. Here SI is a ball in 1" and uq and vq are continuous 
which vanish on the boundary. If uq, vq are nontrivial C 1 functions, the solution 
(u, v) blows up at a finite time T < oo, and ut > 0, vt > on Jl x (0, T), then there 
exist two constants C > c > with 

c . . C c , . C* 

7— 77T- < maXll I,t < 77^ rPTi TT^. wT~ < HiaXll X.t < 7— -77- , 

(T - tfi ~ n (T — ' (T - i)^ " n (T - i)^ ' 

for all i € (0, T). We refer the interested reader, e.g., to [28]. 



2.3 Local lipschitzianity of E 

In the case of a single semilinear elliptic equation, it is known [30] that the ground 
energy map enjoys a basic regularity property, in addition to the continuity, namely 
it is locally Lipschitz continuous (hence differentiable a.e. by virtue of Rademacher's 
theorem). Analogously, for system \S E \ , we obtain the following 

Theorem 2.1 S e Lip loc (R n ). 

Proof. Let p > and (i e R n with \p,\ < p and let be a (dual) solution to (5 M ) 
such that Ifj.^) — (we already know that such a solution does exist, see [1]). 
Then, the corresponding (direct) solution (u M , v^) satisfies 

-Au + u = K(n)v q , -Au + v = Q{p)u p , in R n . (2.7) 

We also know that and are radially symmetric, radially decreasing with 
respect to, say, the origin, and exponentially decaying (see [6,13,27], in particular 
[6, Theorem 2] and [27, Theorem 1(a)]). We claim that there exist vj\ > and 
zu2 = ^U2(po) > independent of /i such that 

roi < ||ip+i < n7 2 and wi < \\Vp\\ L <,+i < w 2 . (2.8) 

Let us prove first the estimates from below. By multiplying the first equation 
of (|2.7|l by u M and taking into account i|1.2|) . we get 

Kill* = / < / 8|K||£ 4+1 ||u / J i , +1 < [3S\\vJ q Lq+ A\^\\m, (2.9) 
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where S is the Sobolev constant. Now, by multiplying the first equation of sys- 
tem (12.711 by and the second equation by and comparing the resulting equa- 
tions, we have 

IKIIU < (f) 9/(9+1 V,ll^ 1)/(9+1) . (2.10) 

By combining inequalities l|2.9|l and (|2.10() , and using again the Sobolev inequality, 
the assertion follows. The proof of the estimate from below for is similar. 

To prove the inequalities from above we simply observe that E is continuous and 

max E(/i) = max I^rf) = max f^u^v^) 
M<po IH<po IH<po 

>(S-5?i)WSSi + (f-p?r)K.ll5ii. 

Thus l|2.8|) follows. As a consequence, according to the definition of the dual norm 
II ' ||^f, we immediately obtain 



aJ wl p + wl q < max |K||^r < P\J ^l P {po) + ^{po)- (2.11) 
v m<po v 



Now, since r? M G A/"^, we get 



(Tif,rp)=l ^4^+/ ^¥^>0. (2.12) 

Hence if S J£+ and, by means of Lemma [2.21 there exists precisely one positive 
number such that Q^ix^rf S By definition, this means that 



(Trfrf) = B( M ,t) 1 ? H-^r / (2.13) 

Moreover, we have 9(/i,fi) = 1. Collecting these facts, we see that, by the implicit 
function theorem, 9 is differentiable with respect to the variable £. Moreover, in 
light of (|2~TT|) . it results that 0(/i,f) remains bounded for /it and £ varying in a 
bounded set. Indeed, by combining (|2.12|l and (|2.13|) . supposing for example that 
p < q, we have 



\m \ p / l%l * 
— / i 
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Then the (local) boundedness of follows immediately by (|2.11|) and by the 



fact that - — - < 0. Let us now observe that 

q p — 



M\ p 



p + 1 7r« qp(c) 



■om^ q 



The gradient of the function {£ i— > 1^(6(^,^)7]^)^ is thus given by 



p-f- 1 



*(m,0 / W,<> 



and so, since 6(]jl,^)7]^ G A/j, in turn we get 
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Kl 



Q— (0 



1^1 



9+1 

A— (0 

Kl p 
Q»(0 



Kl 



■(C) 



9 + 1 



(2.14) 



From this representation formula, the Mean- Value Theorem and the local bound- 
edness of 6, the assertion readily follows. 



2.4 Left and right derivatives of E 

Let us define S(z) as the set of all the positive (dual) solutions (771, 772) of \S z l at the 
energy level £(z)- The representation formulas for the (left and right) directional 
derivatives of £ are provided in the following 
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Theorem 2.2 The directional derivatives from the left and the right ofT, at every 
point z G R" along any w G R™ ea;?s< and it holds 



<9xu 



(z) = sup V z I z (n) ■ w, 



9S W ^ 



inf X7 z I z (n)-w. 

r,£S(z) 



Explicitly, we have 
(*) 



8SV 
<9w / 



sup 

7jes(z) 



<9w 



inf 



1 9Q 

p + 1 <9u> 

1 9Q 



1 dif 



Q^(z) q + ldw 
IjTiI 2 ^ 1 <9if 



(z) 



9+1 

if— (z) 

M ; 

9+1 

K—{z) 



for every z, w 6 R" . 



Proof. Let {/Xj} C R™ be a sequence converging to /xo and let ry- 7 = 7/^ be a 
sequence of (dual) solutions of least energy £(/Zj). We want to prove that, up to a 
subsequence, 

rf* — > 



strongly in Jf, 77 G §(/x ). (2.15) 



Consider the corresponding (direct) solutions (it^ . , . ) (resp. (itoj^o)) of (|2.7[1 
with /x = /ij (resp. /x = (1q). Since (rx^ . , . ) is bounded in VF 2 '(« +1 )/9(R™) x 
VF 2 '( p+1 )/ p (R") (cf. [2]), up to a subsequence, it converges weakly to a pair (uq,vq). 
In addition, since K and Q are uniformly bounded, by virtue of the Schauder local 
regularity estimates (cf. [25]), (u M -,« M ) is bounded in Cj^'f (R") for some (3 > and 



up and Up 



locally in C -sense, 



(2.16) 



so that (u ,vo) solves (|2.7|) with /x = /x . We claim that xto > and «o > 0. 
By [6, Theorem 2], for every j > 1, tx Mj . and are radially symmetric and radially 
decreasing with respect to some point, say the origin, that is 



u N (x)=Uj(r), v H {x) = Vj(r), -^-u 3 -(r) < 0, ^ 3 (r) < 0, 

for every r > 0. Hence, for every j > 1, we have 

u w (0) < -A« w (0) + u H (0) - (0) < [3vl. (0), 

« w (0) < -Av N (0)+v H (0) = Q(m,X,(0) < /3<.(0). 

It follows that, for every j > 1, 

^(0)</3 9+1 <(0). 



(2.17) 
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Then there exists 5 > such that (0) > 6 for every j > 1. Similarly, u^. (0) > <5 
for every j > 1. Hence, letting j — > oo, by (|2.1fi(l . we conclude that uq(0) > 6 
and i>o(0) > <5, which entails uq ^ and i>o ^ 0. Since we have uq > 0, vq > 0, 
K(fi ),Q(fi ) > and 

—Aito + uq > and — Atjo + > 0, 

the claim just follows by a straightforward application of the maximum principle. 
Observe that, by the continuity of E and by Fatou's Lemma, we get 

E( Mo ) = lim E(/*,-) = lim I H (rf) > I^r?) > E(/i ). 

Hence 

lim / w (^) = / M (r ? °) = E( Mo ), 

which reads as 

lim / wj_ = / Um f 9 - f i%i 



j^°°J««-Qp({i j ) JR n Qp(ii ) 3 ^°° J^- n Ki (fij) Jr™ Ki{ijl ) 
In particular, taking into account (|1.2|) . for any <5 > 0, there exists p > such that 



Wi\ p * <s, / i^r« < S, 

{|x|>p} -/{kl>p} 

for every j > 1 sufficiently large. Moreover, of course 



lim / \rji\ p = / 1^1 v , lim / \rf 2 \ « = / |t/ 2 | * . 

3-°° J {\x\<p} J{\x\<p] 3 ^°°J{\x\<p} ->{\x\<P} 

Then we have rfi — > 77° strongly in Jff, namely (|2.15|) holds true. 

Without loss of generality, we can prove the formula of the right derivative of E 
in the case n = 1, z = and w = 1. For any 77 € S(0), we get 

E(p)-E(0)</ p (tf(p,0> ? o )-W) 

= pV e / e (tf(£,0)*7°)|^ e[(M . 

Whence, by virtue of (|2.14|) and the arbitrariness of 77° G §(0), 



hm sup < mf 

P ^o+ P rj°eS(o) 

Moreover, similarly, we get 



Q'(0) /• Iry?!^ tf'(O) f |r$|^ 



P+Ur qV(0) 9 + 1 Jr" if V (0) 



E(p) - E(0) > I p (d{ P , p)tf) - 7 (tf(0, p)^) 
= pV e I e W£,p)r? p )le=/,e[<M, 
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so that, by exploiting i|2.14[) and l|2.15[) . we conclude 

Q'(Q) f Mfg £(o) r m^r - 

P+ 1 Jw Q £ ? L (0) 1 + 1 At« IT^O). ' 

Then the desired formula for the right derivative of E follows. A very similar 
argument provides the corresponding formula for the left derivative. 

Remark 2.1 Nowadays, the further regularity of E is, to our knowledge, an open 
problem. Actually, not even in the case of a single equation the situation is very 
well understood. For instance, on one hand, if we consider the problem 

-e 2 Au + V{x)u = K(x)u p in R" , u>0 in R" , 

then E e C™(R") provided that both the potentials V and K belong to C m (R"), 
with m > 1. On the other hand, if / is not a power (and does not satisfy conditions 
ensuring uniqueness up to translations), for the equation 

~e 2 Au + V{x)u = K{x)f{u) inR™, u>0 inR™, 

we do not know which regularity beyond Lip loc can be achieved by E. Even though 
we do not possess any specific counterexample, our feeling is that there exist func- 
tions / for which the associated E fails to be C 1 smooth. It is evident by the 
(left and right) derivative formulas of E that its further regularity is related to the 
uniqueness of positive radial solutions to — Au + u = f(u) in R™, u > in R™, which 
occurs just for very particular nonlinearities /. Based upon these considerations, 
for semilinear systems, the further regularity of E seems an ever harder matter, 
since as already stressed nothing is known, so far, about the uniqueness of solutions 
to the system 

-Att + u = /(«), -Av + v = g(u), inR™, u, v > inR", 
not even with the particular choices f(v) = v q and g(u) — u p . 

3 Proof of the results 
3.1 Proof of Theorem O 

Let z G £ and let {u £hl v Eh ) £ W 2 '*<r(R n ) x W 2 ' 2 p i (K n ) be a corresponding a 
sequence of strong solutions to \S E l with |u £(i (z)|, |w efc (^)| > S for some S > 0, 
\u £h (z + Ehx)\ — -> 0, \v eh {z + £hx)\ — » as |x| — > oo uniformly w.r.t. h, and 
Sh~ n fe h (u eh ,v eh ) — > E(z) as h ► oo. Let us set: 

<p h (x) = u eh (z + e h x) and i/) h (x) = v £h (z + e h x), 

for all h > 1. Then, since (u £h ,v £h ) is a solution jS £ j , ((phiiph) 1S solution of 

-Aip h + ip h = K(z + e h x)ip q h , -Aiph + i'h = Q{z + £hx)ip p h . (3.1) 



lira inf 

p->0+ 



- g(g) 

P 



> inf 

t7°eS(o) 
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By arguing as in the proof of Theorem 12.21 it is readily proved that, up to a sub- 
sequence, {(p h ) and {ip h ) converge weakly in W 2 ^ q+1 ^ q {M. n ) x W 2 ^ p+1 ^ v {W l ) to 
some (fio and ipo respectively. Let us now prove that there exist 8 > 0, p > and 
ho > 1 such that 

tph(x) < ce- eN and ip h (x) < ce" ^ 1 , for all |x| > p and h > h . (3.2) 

We follow the line of [13]. Since z £ £, then the functions iph and iph decay to zero 
at infinity, uniformly with respect to h. Hence, since p, q > 1, we can find p > 0, 
6 > and ho > 1 such that 

K(z + e h x)^ q h < (1-G 2 )^, 
Q(z + e h x)<p p h < (l-G 2 )^, 

for all |a;| > p and h > ho- Let us set 

!;(x) — pe^ e ^~ p \ p — max max(^ + iph), 

\x\—ph>ho 

and introduce the set 

A = |J D R , 

R>P 

where, for any R > p, we put 

Dr = {p < \x\ < R : iph{x) + ifih(x) > t;(x) for some h > ho}- 
If A = 0, we are done. Instead, if A is nonempty, there exists i?* > p such that 

r P>(n — Lh 

A(e - Vfc - < e 2 - 1 , , j e(») - e 2 ^ - e V 

L \x\ J 

< 6 2 (£-V>h - <fh) < 0, on Dr for all R> R*. 

Hence, by the maximum principle, since — iph ~ fh)\{\x\=p} > 0j we § e t 
% — <Ph > min {0, min (£-iph- <Ph)}, f° r au R> R* 

\x\—R 

so that, letting i? — > oo, yields, for any p > 0, iph(%) + <Ph{x) < f° r I s ! > Pi 
which contradicts the definition of 7^ 0. 

By virtue of the Schauder interior estimates (see e.g. [25]), iph — > 930 and tph — > ipo 
locally in C 2 sense, so that (ipo^o) is a (nontrivial, radial, decaying) solution to 
\S z l - Moreover, in light of the exponential barriers provided by 1)3. 2|) . since z G £, 
it is not difficult to see that (if 0,^0) G ^( z )> f° r we have 



JR» JR™ 

= M<Po,ih) = I*(ti°), 



p+1 
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where rf is the dual solution corresponding to (ipo, ipo). 

Let us now consider the Lagrangian if : ff" x 1 x M x K" x M™ — > M defined as 

££{x, si,s 2 ,Ci,6) = 6 "6 + s i s 2 - 7[^i K ( z + e h x)sl +1 - ^jQ(z + e h x)s\ +1 . 

Then system (|3.1H rewrites as 

-div i^, Vft, V(ph,Viph)) + 9 S2 if (x, <p hl ip h ,V<p h , Vip h ) = 0, in K", 

-div (d^Jif(x,<p h ,iph,\/(ph,Viph)) + d Sl ^(x,(p hl %l> h ,V(p h ,Vij} h ) = 0, in K", 
^, fa > 0, in E™. 

By the Pucci-Serrin identity for systems [21, see §5], we have 

V" / d i q l d^ 2 ) t ^f(x,^ h ,ip hl Vip h ,Vijj h )di'ip h 

n „ 

+ ^2 / d iQ ld (£i)i^{ X 'Ph,i>h,V¥>h,Vi>h)dlLph 

[(div Vh, V^, V^) + q ■ d x 3C(x, <p h ,iph, V^ fc , V-0,0] , 

for all q S C* (R",]R n ). Let us take, for A > 0, 

9 (ar) = (T(Ax),0,...,0), 
and T S C^R") such that T(x) = 1 if |a:| < 1 and T(x) = if \x\ > 2. Then, 

Xd l T(\x)d l if h diiph + V] / \diT(Xx)dii/j h dnph 
= / \diT(\x)^f(x,ip h ,ip h ,Viph,yiph) 

JR™ 

By the arbitrariness of A > 0, letting A — > and keeping h fixed, we obtain 

■^l&iK(z + e h x)ip q h +1 - ^ I d 1 Q(z + = 0. 

Therefore, letting now h — » oo, since in light of (|1.3() we get 

\VK(z + s h x)\,\VQ(z + s h x)\ <ce Meh ^, for |x| large, 
by virtue of (|3.2|) . there holds 



9+ 
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Analogously, we can show that, for all w G 



f \--L-VK(z)^l +1 - ^VQ{z)^ +l 



■w = 0. 



Hence 

^ dw [Z) L (*) 9+1 dw {Z) L K 3 ^ (z) 
Since 77° G §>(z), by Theorem 12. 21 we have 



" Q^(z) yTi ow ' ' ^R- K^r{z) 
Then, by the very definition of (— S)°(z; w) (see Definition I l.ljl . we get 

(-S)°(z; w) > ( ^^^ 0) > 0, for every w G W 1 . 
Then G 9 C (-S)(z) and, since c (-£)(z) = -d c T,{z) (cf. [7]), we obtain z € K. 

3.2 Proof of Corollary fTTTI 

It suffices to combine Theorems 11.11 and 12 .21 taking into account what discussed in 
Section 12 . 21 about the conjectured explicit representation formula for E. 

3.3 Proof of Theorem IPl 

Let m > 1 and z G £ m . The assertion follows by mimicking the various steps in 
the proof of Theorem 11.11 with £ m in place of £ , and combining formula (|3.3|l with 
the definitions of TJ m and /C m , taking into account that rf G Q m (z), as it holds 
Iz(v°) = m i being rf the strong limit of rj £j . Indeed, by Ij3.3|) . there holds 

r+ m (w) > 0, VweM", r- m (tu) > o, Vw e w\ 

so that G ar+ m (0) n 9r- m (0), yielding z G /C m . 
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